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ABSTRACT 

Recent theoretical studies have revealed the possibly important role of the capture and annihilation 
process of weakly interacting massive particles (WIMPs) for the first stars. Using new evolutionary 
models of metal-free massive stars, we investigate the impact of such "dark matter burning" for 
the first stars in different environments of dark matter (DM) halos, in terms of the ambient WIMP 
density (p x ). We find that, in agreement with existing literature, stellar life times can be significantly 
prolonged for a certain range of p x (i.e., 10 10 < p x [GeVcm~ 3 ] < 10 11 with the current upper limit 
for the spin-dependent elastic scattering cross section <Tq D = 5 x 10 _39 cm 2 ). This greatly enhances 
the role of rotationally induced chemical mixing in rotating stars, in favor of abundant production 
of primary nitrogen, massive helium stars and long gamma-ray bursts, from the first stars. We also 
find that stars with p x > 2 x 10 11 GeVcm~ 3 may not undergo nuclear burning stages, confirming the 
previous work, and that ionizing photon fluxes from such DM supported stars are very weak. Delayed 
metal enrichment and slow reionization in the early universe would have resulted if most of the first 
stars had been born in DM halos with such high p x , unless it had been lowered significantly below 
the threshold for efficient DM burning on a short time scale. 

Subject headings: dark matter - early universe - stars: evolution - stars: rotation 



1. INTRODUCTION 

In the current framework of the ACDM cosmologi- 
cal model, the possibility that dark matter (DM) an- 
nihilation effects may impact stellar evolution has re- 
cently received renewed attention. After the origi- 
nal works of Bouquet, Dearborn, Freese, Gould, Gri- 
est, Krauss, Olive, Press, Raffelt, Renzini, Salati, Silk, 
Spergel, Srednicki and Wilczek in the'80s and early '90s, 
several authors have recently re-examined the effects 
that DM, if made of weakly interacting massive particles 
(WIMPs), would have on compact objects (Moskalenko 
& Wai 2007; Bertone & Fairbairn 2008) and on the 
zero-age main sequence of low-mass stars (Fairbairn et 
al. 2008). This exciting activity has been motivated by 
a double scope: any peculiar and distinguishable feature 
of WIMP annihilation on observable stellar quantities is 
extremely precious in the "quest" for dark matter evi- 
dence; on the other hand, all possible effects impacting 
the life of celestial objects must be taken into account by 
astrophysicists in the current precision era. 

In particular, the first stellar episode at high redshift 
occurs under very different conditions from those in the 
present universe. The higher concentration of dark mat- 
ter, the short Hubble time which prevents DM self- 
annihilation from severely affecting the central density, 
and the characteristic formation of a single PopIII star 
in the center of the halo are the most favorable condi- 
tions for DM annihilation effects to be very efficient in 
the first stars. In their pioneering work, Spolyar et al. 



(2008) first found the possibility of very high DM density 
in primordial halos at high redshift to such an extent that 
energy released from DM annihilation at the center may 
halt the gravitational collapse of the baryonic cloud, call- 
ing such a DM powered object a dark star. Iocco (2008, 
hereafter 108) and Freese et al. (2008a) also noticed 
that WIMP capture is most efficient in Population III 
stars. More recently, Freese et al. (2008b, c) and Iocco 
et al. (2008) further investigated the role of annihila- 
tion of adiabatically contracted DM in the formation of 
the first stars. Iocco et al. (2008) have also followed 
the evolution from the pre-main sequence phase to he- 
lium exhaustion in the presence of WIMP capture and 
annihilation (hereafter, DM burning), showing that this 
can severely delay the evolution of pre-MS objects in the 
early universe, as well as extend their MS lifetimes. 

All of these studies motivate us to explore possible 
consequences of DM burning for the final fate of the 
first stars and their feedback effects on the evolution of 
the early universe, even if the role of DM annihilation 
in the formation of the first stars still remains subject 
to many uncertainties (see First Stars III proceedings 
2008, for a recent review). In this Letter, we address 
the issue by discussing the evolution of the first stars 
of 20 < M/Mq < 300 up to the carbon burning stage 
(Sect. 2). We also investigate the interplay of rotation 
with DM burning in the evolution of the first stars of 
100 Mq, given the particular importance of rotation for 
the evolution and deaths of metal poor massive stars (e.g, 
Meynet et al. 2008; Yoon et al. 2008). Implications 
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of our results for the history of reionization in the early 
universe are briefly discussed (Sect. 3). 

2. PHYSICAL ASSUMPTIONS AND RESULTS 

We have implemented the DM capture and annihila- 
tion process in a hydrodynamic stellar evolution code, 
following Gould (1987). The DM capture rate C* is 
calculated using Gould's equations as reported in Eqs. 
(1,2) of 108. Throughout this Letter we assume the DM- 
baryon scattering cross section <jq is 5 x 10~ 39 cm 2 for 
the spin-dependent scattering, to which only hydrogen is 
sensitive, and 10 -43 cm 2 for the spin-independent one. 
These correspond to the current upper limits of WIMP 
direct detection search (Desai et al. 2004; Angle et al. 
2008). We adopt the same values for other parameters 
as those used in 108 to calculate C* , except for the ambi- 
ent DM density (see below). The WIMPs captured by a 
star eventually reach thermal equilibrium with the gas, 
in a configuration dictated by the gravitational poten- 
tial: the consequent DM density can be given by n x {r) — 
n % exp(— r 2 /r 2 ), where r x = c(3kT c /2irp c m x ) 1 / 2 (Griest 
& Seckel 1987). Here T c and p c are temperature and 
density at the stellar center, c and k are the speed of 
light and Boltmann's constant, respectively. The energy 
generation rate due to DM annihilation is given by 



e x (r) = — < crv > n x (r)m x [erg cm 



~ 3 s- 1 ! 



(i) 



where m x is the mass of the DM particle; we assume it 
to be m x = 100 GeV, which is often taken as a fiducial 
value in astrophysical studies for DM search. We use 
< av >= 3 x 10 -26 cm 3 s" 1 , the value best fitting the 
relic DM abundance (see, e.g., Bertone et al. 2005, for a 
recent review). The factor 2/3 is to consider that part of 
the energy is carried away by neutrinos. We consider the 
time dependent evolution of the total number of WIMPs 
(Atot) hi the star, and the number of thermally relaxed 
WIMPs in the core (JVth := / n(r)dV) by the following 
equations in order to normalize n(r): 



dN u 



dt 
dN th 
dt 



C* - J n{r) 2 <av> dV, and (2) 
n(r) 2 < av > dV , (3) 



where r t h is the thermalization rate that can be approx- 
imated by 
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See Iocco et al. (2008) for detailed discussion on the 
thermalization time scale Tth- Our models show that 
equilibrium between C* and Tth is well maintained up to 
the core helium burning phase (see Fig. 2 below). 

We consider several different values for the ambi- 
ent DM energy density p x ranging from to 2 x 
10 12 GeV cm~ 3 . For the initial composition of the 
first stars, we assume the mass fractions of X H, 4 He 
and 3 He to be 0.76, 0.23999, and 0.00001, respectively. 
The mass loss rate from metal free stars is assumed 
to be zero if the Eddington factor Te is smaller than 
0.84 and 10 -14 M yr _1 otherwise, following Krticka & 
Kubat (2006). If a star becomes a helium star by rota- 
tionally induced mixing (see below), we assume that the 
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Fig. 1. — Upper panel: HR diagram of the non-rotating first 
star models on the ZAMS of different masses for different adopted 
values of o"q D p x as indicated by the labels. Here Cq d is the spin 
dependent WIMP scattering cross section, and p x is the ambient 
WIMP density. Stars in the grey shaded region are supposed to be 
only powered by DM burning, without nuclear reactions (i.e., T c < 
10 7 K). We use ct^ d = 5 X 10 -39 cm 2 Lower panel: Life times of the 
non-rotating first stars as a function of the initial mass for different 
values of crg D p x given in the unit of 10 — 26 GeVcm - 1 . These life 
times are obtained from stellar models up to carbon burning for 
Px < 4 x 10l ° GeVcm -3 (i.e.,of D p x < 0.02 X 10 -26 GeVcm -1 ), 
while only approximate estimates are given for p x > 10 11 GeVcm -3 
(i.e.,o-g D p x > 0.05 X 10 -26 GeVcm" 1 ) . 



mass loss rate from such a metal free Wolf-Rayet (WR) 
star is the same as that from a corresponding WR star at 
Z = 10~ 6 , as implied by Vink & de Koeter (2005). The 
code also implements the transport of chemical species 
and angular momentum due to rotationally induced hy- 
drodynamic instabilities and the Spruit-Tayler dynamo. 
Other details on the stellar evolution code are described 
in Yoon et al. (2006). 

Fig. 1 shows the HR diagram of the constructed non- 
rotating first star models on the zero age main sequence 
(ZAMS), for different values of p x . Note that the lumi- 
nosity for a given mass does not significantly change with 
varying p x . This represents the well-known nature of the 
mass-luminosity relation, which is largely independent 
of any particular mode of energy generation (Kippen- 
hahn & Weigert 1990). The stellar structure is adjusted 
such that the ratio of the nuclear to the DM luminosity 
decreases with increasing p x , while the total luminosity 
remains nearly constant, leading to prolonged lifetimes 
of the stars 1 (Fig. 1). If p x is above a critical value 
(px.crit ~ 2 x 10 11 GeVcm" 3 , see Fig. 1) at a given stel- 

1 For a given p x , our models give longer life times than those 
by Iocco et al. (2008). This is because these authors used an 
approximation for the DM capture rate C* , while we integrate Eq 
(2) of 108 over the entire stellar structure. 
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Fig. 2. — Upper Panel: Evolution of the surface luminosity 
(solid line), neutrino luminosity (threc-dotted-dashed line), nuclear 
luminosities due to hydrogen burning (dot-dashed line), helium 
burning (dotted line), carbon burning (long dashed line) and DM 
burning (short dashed line) in the non-rotating 100 Mq model 
sequence. The contribution of DM annihilation to the neutrino 
luminosity, which corresponds to 1/3Ldm, is not included here. 
The rapid increase in Ljj an d Ln c after helium exhaustion is due 
to the hydrogen and helium shell burning. Lower Panel: Evolution 
of the WIMP capture rate (C*, solid line), the thermalization rate 
(r t iu dashed line) and the number of thermalized WIMPs (iVth) 
dotted line) in the corresponding model. See Eqs. (2), (3) and (4). 

lar mass the central density and temperature decrease to 
such an extent that the stars would live forever on the 
ZAMS, having no nuclear reactions. This is in qualitative 
agreement with Iocco et al. (2008) who find that pre-MS 
evolution is "frozen" before reaching the ZAMS at high 
enough p x ; and with Fairbairn et al. (2008), who find 
(for M<4M Q ) that stars move rightward of the ZAMS 
locus on the HR, and eventually join the Hayashi track, 
if "fed" with increasing DM annihilation luminosities. 

No meaningful change in the stellar structure accord- 
ing to different values of p x (< p x , C rit) is observed in non- 
rotating models. Since DM burning only occurs within 
a very small radius r x (<< R CO rc), stars with different 
p x at a given mass have similar amounts of energy flux 
from the core and produce helium cores of a similar size 
(e.g. ~ 40 M from 100 M Q stars). The luminosity 
resulting from DM burning gradually increases early on 
the main sequence as the star expands, but continuously 
decreases in later stages since the significant reduction 
of the number of hydrogen atoms lowers the DM capture 
rate. Rapid increase of the stellar radius after helium ex- 
haustion makes the thermalization time very long, lead- 
ing to reduction of the number of thermalized WIMPs 
(see Fig. 2). The DM luminosity accordingly decreases 
further from 7 x 10 5 L to about 10 5 L Q during the car- 
bon burning phase, in the given example with 100 Mq 
(Fig. 2). Carbon burning and particularly neutrino cool- 
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Fig. 3. — Evolution of the internal structure of rotating 100 Mq 
star models without DM burning (upper panel) and with DM burn- 
ing (lower panel, p x = 4 X 10 10 GeVcm - 3 ). Convective lay- 
ers are hatched, and semi— convective layers are marked by red 
dots. The color shading indicates nuclear energy generation rates. 
The adopted initial rotational velocity at the equatorial surface is 
132 km s — x , which corresponds to 10 % of the Keplerian value. 

ing (L v > 10 10 L Q ) dominate the evolution at this stage 
as shown in Fig 2. As the evolution of the star beyond 
carbon exhaustion should also be governed by neutrino 
cooling and other nuclear reactions such as oxygen burn- 
ing, the effect of DM burning on the pre-supernova struc- 
ture must be minor. As the situation remains similar in 
the other models of 20 < M/M Q < 300, we conclude 
that DM burning may not change the final fate of the 
non-rotating first stars. 

It is noteworthy, however, that rotation can dramati- 
cally change the evolution with DM burning. Fluids in 
rotating massive stars are subject to various rotationally 
induced hydrodynamic instabilities such as Eddingon- 
Sweet circulations, which can cause mixing of chemical 
species across the boundary between the hydrogen burn- 
ing core and the radiative envelope. Such mixing is usu- 
ally stabilized by the strong buoyancy potential due to 
the chemical stratification between the hydrogen burn- 
ing core and the envelope. However, if chemical mix- 
ing can occur faster than the building-up of chemical 
gradients by nuclear burning (i.e., T m i X /r nuc < 1), then 
(quasi-) chemical homogeneity can be maintained on the 
main sequence (so-called chemically homogeneous evolu- 
tion) (Maeder 1987). The condition, T mix /T nuc < 1, can 
be met either by reducing r m i x or increasing r nuc . Rapid 
rotation tends to do the former, and we find that DM 
burning tends to do the latter. This is because the mixing 
time scale due to Eddington-Sweet circulations remains 
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almost unchanged with increasing p x (for p x < p x . cr it), 
while the nuclear burning time scale increases due to the 
DM burning. This explains the remarkable impact of DM 
burning in the evolution of our rotating models as shown 
in Fig. 3. As shown in the upper panel, rotation in this 
model (winit/wKopior = 0.1) is not rapid enough to cause 
the strong mixing by itself without DM burning. With 
p x = 4 x 10 10 GeV cm -3 (see lower panel), however, the 
100 M Q star lives about 10 times longer than in the corre- 
sponding non-DM burning case, and the star undergoes 
the quasi- chemically homogeneous evolution even with 
slow initial rotation because of the DM burning effect on 
nuclear burning time scale (cf. Yoon et al. 2006). The 
star is thus gradually transformed into a massive helium 
star by the end of main sequence as almost all of the 
hydrogen atoms in the star are fused into heliums due to 
mixing. 

3. DISCUSSION 

Our results indicate that DM burning should not sig- 
nificantly alter our view on the final fate of the non 
rotating first stars: pair-instability supernovae for 140 < 
M/Mq < 260, and core-collapse events for other masses 
(Heger & Woosley 2002), although their life times may 
be significantly prolonged. However, the impact of DM 
burning appears more important for rotating stars. The 
quasi-chemically homogeneous evolution can be rather 
easily realized even with moderate rotation velocities 
with 10 10 < p x [GeVcm- 3 ] < 10 11 . Such evolution can 
lead to production of massive helium stars that emit large 
amounts of helium ionizing photons, as well as abundant 
production of primary nitrogen, as discussed in Yoon et 
al. (2008). Note also that the quasi-chemically homo- 
geneous evolution scenario (CHES) is one of the favored 
ones for the production of long GRBs from metal poor 
stars (Yoon & Langer 2005; Woosley & Heger 2006). 
Our result therefore indicates that DM burning might 
promote the production of long gamma-ray bursts from 
the first stars of 12 < M/M < 60 via the CHES channel 
(see Yoon et al. 2006). 

DM burning in the first stars must have consequences 
in the history of reionization in the early universe. Ta- 
ble 1 lists the number of hydrogen and helium ioniz- 
ing photons emitted from 100 M models. If p x < 
2 x 10 11 GeVcm -3 , the total number of ionizing photons 
increases proportionally to the DM density, as a direct 
consequence of the life-prolonging effect of DM burning. 
For a given p x rotation does not significantly alter hy- 
drogen ionizing photon counts nor the lifetime. However, 
for models that take path of the CHES, the helium ion- 
izing photon counts is increased by more than a factor 
of 2 compared to the non-rotating case. If p x is very 
large, on the other hand, the surface temperature of the 
star drops significantly enough (Fig. 1) that the total 
number of ionizing photons is reduced even with much 
longer lifetimes due to the DM burning. For example, 
with p x = 2 x 10 12 GeVcm" 3 , it would have to take ~ 
10 Gyr to emit as many hydrogen ionizing photons as a 
non-DM burning counterpart. Therefore, if most of the 
first stars had been born with such high p x their con- 
tribution to reionization would have been dramatically 



reduced. The effect of the temperature drop on the num- 
ber of helium ionizing photons is even more prominent: 
with p x = 2 x 10 12 GeVcm -3 , it decreases by more than 

TABLE 1 

Number of hydrogen and helium ionizing photons from 
100 m star models 



v^tpUK Px [GeVcm 3 ] TVjj JVh c Duration 



0.0 O00 1.2 X 10 M 2.2 X 10 b2 3.2 Myr 

0.0 2 x 10 10 2.1 x 10 64 3.4 x 10 62 5.5 Myr 

0.0 4 X 10 10 8.5 X 10 64 1.5 X 10 63 22.3 Myr 

0.0 10 11 2.9 X 10 65 6.3 X 10 62 100.0 Myr* 

0.0 2 X 10 11 2.0 X 10 65 1.7 X 10 61 100.0 Myr* 

0.0 2 X 10 12 1.3 X 10 62 9.4 X 10 43 100.0 Myr* 

Ol O00 1.5 X 10 04 2.5 X 10° 2 3.4 Myr 

0.1 2 x 10 10 2.7 x 10 64 5.2 x 10 62 6.0 Myr 

0.1 4 X 10 1(1 8.7 X 10 64 3.8 X 10 63 19.6 Myr 



The numbers are calculated only for the first 100 Myr. 

19 orders of magnitude compared to the other cases with 
lower p x . Future study on the history of helium ioniza- 
tion at high redshift might therefore be a strong probe 
of DM burning in the first stars. 

In this study we assume that the background DM 
halo density stays constant throughout the stellar evo- 
lution. This assumption may be valid if the stellar life 
times are shorter than about 100 Myr - which is the ex- 
pected merger timescale of DM halos at at z <~ 20 (e.g., 
Lacey & Cole 1993) - as in our model sequences with 
Px ~ 4x 10 10 GeVcm -3 . The evolution of DM burners in 
halos with higher p x , however, should be critically deter- 
mined by the change of DM halo environments. If p x is 
sufficiently reduced due to merger events and/or to dis- 
placement of the star from the densest region of the DM 
halo, the DM burners will become "normal" stars, domi- 
nated by nuclear burning. The star may then die quickly 
as we expect for normal stars, or become a "born-again" 
DM burner if p x increases again in later stages for some 
reason. The detailed history of the feedback from the 
first stars (e.g. metal enrichment and reionization) on 
the evolution of the early universe may depend on the 
nature and evolution of DM halos where the first stars 
are formed. This issue should be addressed in future 
work. 

As a note added in proof, we acknowledge that Taoso et 
al. (2008) independently report similar results about the 
effect of DM burning on the MS life time and some stellar 
properties, using a different numerical code. S. A. appre- 
ciates helpful discussion with M. Alvarez and M.Busha. 
S. A. is supported by a Department of Energy Contract 
to SLAC DE-AC3-76SF00515. F. I. is grateful to P. 
Scott and M. Taoso for helpful discussion. F. I. is sup- 
ported by MIUR through grant PRIN-2006. S.C. Y. 
is grateful to C. Church for the help with the text, 
and to W. Hillebrandt and S. Woosley for supporting 
his visit to MPA, Garching, in June, 2008, where part 
of the manuscript has been prepared. S.C. Y. is sup- 
ported by the DOE SciDAC Program (DOE DE-FC02- 
06ER41438). 
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